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Abstract. 

I review recent observational and theoretical progress in our understanding of the 
cosmic evolution of luminous sources. Largely due to a combination of deep HST imag- 
ing, Keck spectroscopy, and COBE far-IR background measurements, new constraints 
have emerged on the emission history of the galaxy population as a whole. Barring 
large systematic effects, the global ultraviolet, optical, near- and far-IR photometric 
properties of galaxies as a function of cosmic time cannot be reproduced by a simple 
stellar evolution model defined by a constant (comoving) star-formation density and 
a universal (Salpeter) IMF, and require instead a substantial increase in the stellar 
birthrate with lookback time. While the bulk of the stars present today appears to 
have formed relatively recently, the existence of a decline in the star-formation density 
above z ~ 2 remains uncertain. The history of the transition from the cosmic 'dark 
age' to a ionized universe populated with luminous sources can constrain the star for- 
mation activity at high redshifts. If stellar sources are responsible for photoionizing the 
intergalactic medium at z as 5, the rate of star formation at this epoch must be com- 
parable or greater than the one inferred from optical observations of galaxies at z ~ 3. 
A population of dusty, Type II AGNs at z < 2 could make a significant contribution 
to the FIR background if the accretion efficiency is ~ 10%. 



INTRODUCTION 

There is little doubt that the last few years have been exciting times in galaxy 
formation and evolution studies. The remarkable progress in our understanding of 
faint galaxy data made possible by the combination of HST deep imaging [62] and 
ground-based spectroscopy [35], [14], [56] has permitted to shed new light on the 
evolution of the stellar birthrate in the universe, to identify the epoch 1 < z < 2 
where most of the optical extragalactic background light was produced, and to 
set important contraints on galaxy evolution scenarios [41], [57], [3], [23]. The ex- 
plosion in the quantity of information available on the high-redshift universe at 
optical wavelengths has been complemented by the detection of the far-IR/sub- 
mm background by DIRBE and FIRAS [28], [16]. The IR data have revealed the 
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FIGURE 1. Spectrum of the extragalactic background light as derived from a compilation of 
ground-based and space-based galaxy counts in the U, J3, V, I, and if-bands (filled dots), together 
with the FIRAS 125-5000 /im (solid and dashed lines) and DIRBE 140 and 240 /im (filled squares) 
detections. The empty squares show the DIRBE points after correction for WIM dust emission 
[32]. 



optically 'hidden' side of galaxy formation, and shown that a significant fraction of 
the energy released by stellar nucleosynthesis is re-emitted as thermal radiation by 
dust. The underlying goal of all these efforts is to understand the growth of cosmic 
structures and the mechanisms that shaped the Hubble sequence, and ultimately 
to map the transition from the cosmic 'dark age' to a ionized universe populated 
with luminous sources. While one of the important questions recently emerged is 
the nature (starbursts or AGNs?) and redshift distribution of the ultraluminous 
sub- mm sources discovered by SCUBA [30], [2], [36], of perhaps equal interest is the 
possible existence of a large population of faint galaxies still undetected at high-z, 
as the color-selected ground-based and Hubble Deep Field (HDF) samples include 
only the brightest and bluest star-forming objects. In hierarchical clustering cos- 
mogonies, high- z dwarfs and/or mini-quasars (i.e. an early generation of stars and 
accreting black holes in dark matter halos with circular velocities v c ~ 50 kms" 1 ) 
may actually be one of the main source of UV photons and heavy elements at early 
epochs [44], [25], [26]. 

In this talk I will focus on some of the open issues and controversies surround- 
ing our present understanding of the history of the conversion of cold gas into 
stars within galaxies, and of the evolution with cosmic time of luminous sources 
in the universe. An Einstein-deSitter (EdS) universe (Qm — 1 ; = 0) with 
h = H /100 kms -1 Mpc -1 = 0.5 will be adopted in the following. 



OPTICAL/FIR BACKGROUND 



The extragalactic background light (EBL) is an indicator of the total luminosity 
of the universe. It provides unique information on the evolution of cosmic struc- 
tures at all epochs, as the cumulative emission from galactic systems and AGNs is 
expected to be recorded in this background. Figure 1 shows the optical EBL from 
known galaxies together with the recent COBE results. The value derived by inte- 
grating the galaxy counts [49] down to very faint magnitude levels [because of the 
flattening at faint magnitudes of the N(m) differential counts most of the contribu- 
tion to the optical EBL comes from relatively bright galaxies] implies a lower limit 
to the EBL intensity in the 0.3-2.2 fxm interval of J opt ~ 12 nWm~ 2 sr -1 . 1 When 
combined with the FIRAS and DIRBE measurements (/fir ~ 16 nWm^sr 1 
in the 125-5000 /iin range), this gives an observed EBL intensity in excess of 
28 nWm~ 2 sr _1 . The correction factor needed to account for the residual emis- 
sion in the 2.2 to 125 jum region is probably < 2 [13]. We shall see below how a 
population of dusty AGNs could make a significant contribution to the FIR back- 
ground. In the rest of this talk I will adopt a conservative reference value for the 
total EBL intensity associated with star formation activity over the entire history 
of the universe of Jebl = 40 ho nWm -2 sr _1 . 

COSMIC STAR FORMATION 

It has become familiar to interpret recent observations of high-redshift sources 
via the comoving volume-averaged history of star formation. This is the mean 
over cosmic time of the stochastic, possibly short-lived star formation episodes of 
individual galaxies, and follows a relatively simple dependence on redshift. Its 
latest version, uncorrected for dust extinction, is plotted in Figure 2 (left). The 
measurements are based upon the rest-frame UV luminosity function (at 1500 and 
2800 A), assumed to be from young stellar populations [38]. The prescription for 
a 'correct' de- reddening of these values has been the subject of an ongoing debate. 
Dust may play a role in obscuring the UV continuum of Canada-France Reshift 
Survey (CFRS, 0.3 < z < 1) and Lyman-break (z ~ 3) galaxies, as their colors are 
too red to be fitted with an evolving stellar population and a Salpeter initial mass 
function (IMF) [41]. The fiducial model of [41] had an upward correction factor 
of 1.4 at 2800 A, and 2.1 at 1500 A. Much larger corrections have been argued 
for by [50] (xlO at z = 1), [43] (xl5 at z = 3), and [51] (xl6 at z > 2). As 
noted already by [38] and [41], a consequence of such large extinction values is the 
possible overproduction of metals and red light at low redshifts. Most recently, 
the evidence for more moderate extinction corrections has included measurements 
of star-formation rates (SFR) from Balmer lines by [59] (x2 at z = 0.2), [21] 

Note that the direct detection of the optical EBL at 3000, 5500, and 8000 A derived from HST 
data by [5] implies values that are about a factor of two higher than the integrated light from 
galaxy counts. 
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FIGURE 2. Left: Mean comoving density of star formation as a function of cosmic time. The 
data points with error bars have been inferred from the UV-continuum luminosity densities of [35] 
(filled dots), [10] (filled squares), [41] (filled pentagons), [60] (empty dot), and [58] (empty square). 
The dotted line shows the fiducial rate, (p\) = 0.054 M Q yr _1 Mpc -3 , required to generate the 
observed EBL. Right: dust corrected values (A 150 o = 1.2 mag, SMC-type dust in a foreground 
screen). The Ha determinations of [19], [59], and [21] (filled triangles), together with the SCUBA 
lower limit [30] (empty pentagon) have been added for comparison. 



(x3.1 ± 0.4 at z = 1), and [48] (x2-6atz = 3). ISO follow-up of CFRS fields 
[17] has shown that the star-formation density derived by FIR fluxes (x2.3 ± 0.7 
at < z < 1) is about 3.5 times lower than in [50]. Figure 2 (right) depicts an 
extinction-corrected (with y4 150 o = 1.2 mag, 0.4 mag higher than in [41]) version of 
the same plot. The best-fit cosmic star formation history (shown by the dashed- 
line) with such a universal correction produces a total EBL of 37 nWm _2 sr _1 . 
About 65% of this is radiated in the UV+optical+near-IR between 0.1 and 5 /j,m; 
the total amount of starlight that is absorbed by dust and reprocessed in the far-IR 
is 13 nWm -2 sr _1 . Because of the uncertainties associated with the incompleteness 
of the data sets, photometric redshift technique, dust reddening, and UV-to-SFR 
conversion, these numbers are only meant to be indicative. On the other hand, this 
very simple model is not in obvious disagreement with any of the observations, and 
is able, in particular, to provide a reasonable estimate of the galaxy optical and 
near-IR luminosity density. 

STELLAR BARYON BUDGET 

With the help of some simple stellar population synthesis tools it is possible 
at this stage to make an estimate of the stellar mass density that produced the 
integrated light observed today. The total bolometric luminosity of a simple stel- 
lar population (a single generation of coeval stars) having mass M can be well 
approximated by a power-law with time for all ages t > 100 Myr, 
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(cf. [8]), where we have assumed solar metallicity and a Salpeter IMF truncated 
at 0.1 and 125 M . In a stellar system with arbitrary star-formation rate per unit 
cosmological volume, p*, the comoving bolometric emissivity at time t is given by 
the convolution integral 

Pboi(t) = [* L(r)p*(t-T)dT. (2) 
Jo 

The total background light observed at Earth (t = £#) is 

T C ft* p bo i(t) 

Jebl = — / — — at, (3) 

47T 7o 1 + z 

where the factor (1 + z) at the denominator is lost to cosmic expansion when 
converting from observed to radiated (comoving) luminosity density. From the 
above equations it is easy to derive in a EdS cosmology 

7 EBL = 740 n Wm -V-( MsyT f; Mpc _, > (^) ■ (4) 

The observations shown in Figure 1 therefore imply a "fiducial" mean star for- 
mation density of (p*) = 0.054 J 40 MQyr" 1 Mpc~ 3 . In the instantaneous recycling 
approximation, the total stellar mass density observed today is 

P*(t H ) = (1 - R) f tH p*(t)dt w 5 x 10 8 1 40 M Mpc~ 3 (5) 
Jo 

(corresponding to = 0.007 7 40 ), where R is the mass fraction of a generation of 
stars that is returned to the interstellar medium, R ks 0.3 for a Salpeter IMF. The 
optical/FIR background therefore requires that about 10% of the nucleosynthetic 
baryons today [7] are in the forms of stars and their remnants. The predicted 
stellar mass-to-blue light ratio is (M/Lb) ~ 5. These values are quite sensitive to 
the lower-mass cutoff of the IMF, as very-low mass stars can contribute significantly 
to the mass but not to the integrated light of the whole stellar population. A lower 
cutoff of 0.5 M Q instead of the 0.1 M Q adopted would decrease the mass-to-light 
ratio (and fi*) by a factor of 1.9 for a Salpeter function. 



TWO SIMPLE MODELS 



Based on the agreement between the z ~ 3 and z ~ 4 luminosity functions at 
the bright end, it has been recently argued [58] that the decline in the luminosity 
density of faint HDF Lyman-break galaxies observed in the same redshift interval 
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FIGURE 3. Left: Synthetic evolution of the near-IR luminosity density at rest-frame wave- 
lengths of 1.0 (long-dashed line) and 2.2 fim (short-dashed line). The model assumes a constant 
star-formation rate of p* = 0.054 M yr" 1 Mpc~ 3 (Salpeter IMF). The dotted (2.2 fxm) and 
dash- dotted (1.0 /jm) curves show the emissivity of a simple stellar population with formation 
redshift z on = 5, and total mass equal to the mass observed in spheroids today [18]. The data 
points are taken from [35] (filled dots) and [20] (filled triangle). Right: Same but in the B-band. 
The data points are taken from [35] (filled dots), [14] (empty triangles), and [10] (empty squares). 

[38] may not be real, but simply due to sample variance in the HDF. When ex- 
tinction corrections are applied, the emissivity per unit comoving volume due to 
star formation may then remain essentially flat for all redshift z > 1 (see Fig. 2). 
While this has obvious implications for hierarchical models of structure formation, 
the epoch of first light, and the reionization of the intergalactic medium (IGM), 
it is also interesting to speculate on the possibility of a constant star-formation 
density at all epochs < z < 5, as recently advocated by [47]. Figure 3 shows 
the time evolution of the blue and near-IR rest-frame luminosity density of a stel- 
lar population characterized by a Salpeter IMF, solar metallicity, and a (constant) 
star- formation rate of p* = 0.054 M yr -1 Mpc -3 (needed to produce the observed 
EBL). The predicted evolution appears to be a poor match to the observations: it 
overpredicts the local B and .fT-band luminosity densities, and underpredicts the 
1 /im emissivity at z ~ 1 from the CFRS survey. 2 

At the other extreme, we know from stellar population studies that about half 
of the present-day stars are contained in spheroidal systems, i.e. elliptical galaxies 
and spiral galaxy bulges, and that these stars formed early and rapidly [4]. The 
expected rest-frame blue and near-IR emissivity of a simple stellar population with 
formation redshift z on = 5 and total mass density equal to the mass in spheroids 

2 ) The near-IR light is dominated by near-solar mass evolved stars, the progenitors of which make 
up the bulk of a galaxy's stellar mass, and is more sensitive to the past star-formation history 
than the blue light. 



observed today (see below) is shown in Figure 3. HST-NICMOS deep observations 
may be able to test similar scenarios for the formation of elliptical galaxies at early 
times. 

TYPE II AGNS 

Recent dynamical evidence indicates that supermassive black holes reside at the 
center of most nearby galaxies. The available data (about 30 objects) show a strong 
correlation (but with a large scatter) between bulge and black hole mass [42] , with 
= 0.006 Mbuigc as a best-fit. The total mass density in spheroids today is 
f^buige = 0.0036to;ooi7 [18], implying a mean mass density of dead quasars 

p bh = 1.5^;° 7 xlO 6 M Mpc- 3 . (6) 

Noting that the observed energy density from all quasars is equal to the emitted 
energy divided by the average quasar redshift [63], the total contribution to the 
EBL from accretion onto black holes is 

4h = j-tt^t ~ 18 nWm- 2 sr-Vi(l + z)~\ (7) 

where 770.1 is the efficiency for transforming accreted rest-mass energy into radia- 
tion (in units of 10%). A population of AGNs at (say) z ~ 1.5 could then make 
a significant contribution to the FIR background if dust-obscured accretion onto 
supermassive black holes is an efficient process [24], [15]. It is interesting to note 
in this context that a population of AGNs with strong intrinsic absorption (Type 
II quasars) is actually invoked in many current models for the X-ray background 
[37], [9]. 

SOURCES OF IONIZING RADIATION 

The application of the Gunn-Peterson constraint on the amount of smoothly 
distributed neutral material along the line of sight to distant objects requires the 
hydrogen component of the diffuse IGM to have been highly ionized by z ~ 5 
[53], and the helium component by z ~ 2.5 [12]. From QSO absorption studies we 
also know that neutral hydrogen at early epochs accounts for only a small fraction, 
~ 10%, of the nucleosynthetic baryons [33]. It thus appears that substantial sources 
of ultraviolet photons were present at z <: 5, perhaps low-luminosity quasars [26] 
or a first generation of stars in virialized dark matter halos with T V - 1T ~ 10 4 — 
10 5 K [46], [25], [44]. Early star formation provides a possible explanation for the 
widespread existence of heavy elements in the IGM [11], while reionization by QSOs 
may produce a detectable signal in the radio extragalactic background at meter 
wavelengths [40]. Establishing the character of cosmological ionizing sources is an 




redshift redshift 

FIGURE 4. Left: comoving space density of bright QSOs as a function of redshift. The data 
points with error bars are taken from [27] (filled dots), [61] (filled squares), [52] (crosses), and 
[31] (filled pentagon). The empty triangles show the space density of the Parkes flat-spectrum 
radio- loud quasars with P > 7.2 x 10 26 W Hz -1 sr _1 [29]. Right: comoving emission rate of 
hydrogen Lyman-continuum photons (solid line) from QSOs, compared with the minimum rate 
(dashed line) which is needed to fully ionize a fast recombining (with ionized gas clumping factor 
C = 30) EdS universe with £lbh 2 = 0.02. Models based on photoionization by quasar sources 
appear to fall short at z ~ 5. The data point shows the estimated contribution of star-forming 
galaxies at z *» 3, assuming that the fraction of Lyman continuum photons which escapes the 
galaxy H i layers into the intergalactic medium is / csc = 0.5 (see [39] for details). 



efficient way to constrain competing models for structure formation in the universe, 
and to study the collapse and cooling of small mass objects at early epochs. 

What keeps the universe ionized at z = 5? The problem can be simplified by 
noting that the breakthrough epoch (when all radiation sources can see each other 
in the Lyman continuum) occurs much later in the universe than the overlap epoch 
(when individual ionized zones become simply connected and every point in space 
is exposed to ionizing radiation). This implies that at high redshifts the ionization 
equilibrium is actually determined by the instantaneous UV production rate [39]. 
The fact that the IGM is rather clumpy and still optically thick at overlapping, 
coupled to recent observations of a rapid decline in the space density of radio-loud 
quasars and of a large population of star- forming galaxies ai z > 3, has some 
interesting implications for rival ionization scenarios and for the star formation 
activity at < 3 < z < 5. 

The existence of a decline in the space density of bright quasars at redshifts 
beyond ~ 3 was first suggested by [45], and has been since then the subject of 
a long-standing debate. In recent years, several optical surveys have consistently 
provided new evidence for a turnover in the QSO counts [27], [61], [52], [31]. The 
interpretation of the drop-off observed in optically selected samples is equivocal, 
however, because of the possible bias introduced by dust obscuration arising from 
intervening systems. Radio emission, on the other hand, is unaffected by dust, and 
it has recently been shown [54] that the space density of radio-loud quasars also 
decreases strongly for z > 3. This argues that the turnover is indeed real and that 
dust along the line of sight has a minimal effect on optically-selected QSOs (Figure 
4, left). The QSO emission rate (corrected for incompleteness) of hydrogen ionizing 
photons per unit comoving volume is shown in Figure 4 (right) [39]. 

Galaxies with ongoing star-formation are another obvious source of Lyman con- 
tinuum photons. Since the rest-frame UV continuum at 1500 A (redshifted into the 
visible band for a source at z ~ 3) is dominated by the same short-lived, massive 
stars which are responsible for the emission of photons shortward of the Lyman 
edge, the needed conversion factor, about one ionizing photon every 10 photons at 
1500 A, is fairly insensitive to the assumed IMF and is independent of the galaxy 
history for t ^> 10 7 yr. Figure 4 (right) shows the estimated Lyman-continuum lu- 
minosity density of galaxies at z ~ 3. 3 The data point assumes a value of f esc = 0.5 
for the unknown fraction of ionizing photons which escapes the galaxy H 1 layers 
into the intergalactic medium. A substantial population of dwarf galaxies below 
the detection threshold, i.e. having star-formation rates < 0.3 M yr _1 , and with a 
space density in excess of that predicted by extrapolating to faint magnitudes the 
best-fit Schechter function, may be expected to form at early times in hierarchical 
clustering models, and has been recently proposed by [44] and [39] as a possible 
candidate for photoionizing the IGM at these epochs. One should note that, while 
highly reddened galaxies at high redshifts would be missed by the dropout color 



3 ) At all ages > 0.1 Gyr one has L(1500)/L(912) w 6 for a Salpeter mass function and constant 
SFR [6]. This number neglects any correction for intrinsic H 1 absorption. 



technique (which isolates sources that have blue colors in the optical and a sharp 
drop in the rest-frame UV), it seems unlikely that very dusty objects (with / esc <C 1) 
would contribute in any significant manner to the ionizing metagalactic flux. 



REIONIZATION OF THE IGM 

When an isolated point source of ionizing radiation turns on, the ionized volume 
initially grows in size at a rate fixed by the emission of UV photons, and an ion- 
ization front separating the H n and H i regions propagates into the neutral gas. 
Most photons travel freely in the ionized bubble, and are absorbed in a transition 
layer. The evolution of an expanding H n region is governed by the equation 

dV - SHV^^-YL, (8) 



dt 71h t 



rcc 



where V} is the proper volume of the ionized zone, iVi on is the number of ionizing 
photons emitted by the central source per unit time, % is the mean hydrogen 
density of the expanding IGM, H is the Hubble constant, and t rcc is the hydrogen 
mean recombination timescale, 

t rec = [(l + 2 X )n H aBC]- 1 = 0.3Gyr^) '(^fVo 1 - (9) 

One should point out that the use of a volume-averaged clumping factor, C, in the 
recombination timescale is only justified when the size of the H n region is large 
compared to the scale of the clumping, so that the effect of many clumps (filaments) 
within the ionized volume can be averaged over (see Figure 5). Across the I-front 
the degree of ionization changes sharply on a distance of the order of the mean free 
path of an ionizing photon. When t rcc <C t, the growth of the H n region is slowed 
down by recombinations in the highly inhomogeneous medium, and its evolution 
can be decoupled from the expansion of the universe. Just like in the static case, the 
ionized bubble will fill its time- varying Stromgren sphere after a few recombination 
timescales, 

V I = ] ^^(l-e- t ^). (10) 
n H 

In analogy with the individual H n region case, it can be shown that hydrogen 
component in a highly inhomogeneous universe is completely reionized when the 
number of photons emitted above 1 ryd in one recombination time equals the mean 
number of hydrogen atoms [39] . At any given epoch there is a critical value for the 
photon emission rate per unit cosmological comoving volume, 
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FIGURE 5. Propagation of an ionization front in a 128 3 cosmological density field produced 
by a mini-quasar with N — 5 x 10 53 s -1 . The box length is 2.4 comoving Mpc. The quasar is 
turned on at the densest cell, which is found in a virialized halo of total mass 1.3 x 10 11 Mq. The 
solid contours give the position of the I-front at 0.15, 0.25, 0.38, and 0.57 Myr after the quasar 
has switched on at z = 7. The underlying greyscale image indicates the initial H i density field. 
(From [1].) 

independently of the (unknown) previous emission history of the universe: only 
rates above this value will provide enough UV photons to ionize the IGM by that 
epoch. One can then compare our estimate of Afi on to the estimated contribution 
from QSOs and star-forming galaxies. The uncertainty on this critical rate is 
difficult to estimate, as it depends on the dumpiness of the IGM (scaled in the 
expression above to the value inferred at z = 5 from numerical simulations [22]) 
and the nucleosynthesis constrained baryon density The evolution of the critical 
rate as a function of redshift is plotted in Figure 4 (right). While A/ion is comparable 
to the quasar contribution at z > 3, there is some indication of a deficit of Lyman 
continuum photons at z — 5. For bright, massive galaxies to produce enough UV 
radiation at z — 5, their space density would have to be comparable to the one 
observed at z ~ 3, with most ionizing photons being able to escape freely from 
the regions of star formation into the IGM. This scenario may be in conflict with 
direct observations of local starbursts below the Lyman limit showing that at most 
a few percent of the stellar ionizing radiation produced by these luminous sources 
actually escapes into the IGM [34]. 4 

It is interesting to convert the derived value of A/i on into a "minimum" SFR per 
unit (comoving) volume, p* (hereafter we assume VL^h 2 = 0.02 and C = 30): 

p*(z)=Mon(z) x lO- 53 ' 1 ^ 1 « 0.013/^ (^) 3 M yr- 1 Mpc- 3 . (12) 

4 ) Note that, at z = 3, Lyman-break galaxies would radiate more ionizing photons than QSOs 
for / csc > 30%. 



The star formation density given in the equation above is comparable with the 
value directly "observed" (i.e., uncorrected for dust reddening) at z 3 [41]. The 
conversion factor assumes a Salpeter IMF with solar metallicity, and has been 
computed using a population synthesis code [6]. It can be understood by noting 
that, for each 1 M & of stars formed, 8% goes into massive stars with M > 2OM 
that dominate the Lyman continuum luminosity of a stellar population. At the end 
of the C-burning phase, roughly half of the initial mass is converted into helium 
and carbon, with a mass fraction released as radiation of 0.007. About 25% of the 
energy radiated away goes into ionizing photons of mean energy 20 eV. For each 1 
M Q of stars formed every year, we then expect 



0.08 x 0.5 x 0.007 x 0.25 x M c 2 1 
20 eV Tyr 

to be emitted shortward of 1 ryd. 



10 53 phots" 1 (13) 
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